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Accretion disk onto boson stars: a way to supplant black holes candidates 
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The emission spectrum from a simple accretion disk model around a compact object is compared 
for the cases of a black hole (BH) and a boson star (BS) playing the role of the central object. 

It was found in the past that such a spectrum presents a hardening at high frequencies; however, 
here it is shown that the self-interaction and compactness of the BS have the effect of softening 
the spectrum, the less compact the star is, the softer the emission spectrum at high frequencies. 
Because the mass of the boson fixes the mass of the star and the self-interaction the compactness 
of the star, we find that, for certain values of the BS parameters, it is possible to produce similar 
spectra to those generated when the central object is a BH. This result presents two important 
implications: (i) using this simple accretion model, a BS can supplant a BH in the role of compact 
object accreting matter, and (ii) within the assumptions of the present accretion disk model we do 
not find a prediction that could help distinguish a BH from a BS with appropriate parameters of 
mass and self-interaction. 

PACS numbers: 97.10.Gz 04.40.-b 04.70.-s 98.35. Jk 


Since boson stars (BSs) were found to be stationary 
solutions of the Einstein-Klein-Gordon (EKG) system of 
equations Hi there have been two approaches for their 
study in astrophysical grounds: i) proof their existence 
through measurable predictions and ii) use BSs as toy 
models for problems requiring an exotic matter compo¬ 
nent. Within the first aim some topics have been devel¬ 
oped: prediction of particular signatures in gravitational 
wave emission of perturbed boson stars illii ; pre¬ 
diction of differences in the emission spectrum in accre¬ 
tion disks around supermassive BSs vs that from BHs 
Hi! The second line of research involves propos¬ 
ing gigantic BSs as galactic dark halos 0, (Til, US! 
BSs as alternative to MACHOS 14] and BSs as solu¬ 
tions of tensor-scalar theories mnim Around these 
two fronts the understanding of BSs has been widely de¬ 
veloped, and the most important results can be found in 
two recent reviews by Jetzer 0 and Schunck and Mielke 
& In this manuscript we focus on the first approach, 
and try to provide a picture of the properties that make 
a BS look like a BH, which we consider to be important 
because observations are at the verge of proving -or not- 
that black hole candidates (BHCs) have a horizon or not 
and when such debate involves BSs ill- Firstly 
we present the fundamentals on Boson Stars as follows. 

Boson stars are stationary solutions of Einstein’s equa¬ 
tions G^ = 8-kGT A, where = ±[S M </>* d v (/> + 

d^dvcj)*} - + V(\4>\ 2 )\ is the stress-energy 

tensor of a complex scalar field (j>, g^ is the metric 
of the space-time and V is the potential of the field, 
which for the present case is restricted to have the form 
U(|0|) = ^m 2 \(j)\ 2 + j\(/)\ 4 where m is the mass of the field 
and A its self-interaction. Whether or not the scalar field 
plays the role of a fundamental spinless particle, one can 
always consider it appears in an effective theory with la- 
grangian density C = — 16 ^ G +g Mi 'd /x 0*d t ,</>+ 1 F(|</>|). This 
mean field approach helps at avoiding the formulation of 


quantum field theory in a curved background. 

Although axi-symmetric BS solutions have been found, 
we restrict here to the spherically symmetric case and 
choose the line element of space-time in Schwarzscliild 
coordinates to be ds 2 = —a{r) 2 dt 2 + a(r) 2 dr 2 + r 2 dfl 2 . 
Assuming the field to have a time dependence <p(r, t) = 
(j)o{r)e~ lut the stress energy tensor is time independent, 
which implies that the space-time is stationary and there¬ 
fore that the metric functions a and a depend only on the 
radial coordinate r. Under this assumptions the Einstein- 
Klein-Gordon equations are: 
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where we have used the rescaled variables <j >o —> V4ttG(/) 0 , 
r — > mr, t — > tot, a —> (j a and A = A 7 T Q m 2 ; primes indi¬ 
cate derivative with respect to the radial coordinate. The 
system (1-3) is a set of coupled ordinary differential equa¬ 
tions to be solved under the conditions a(0) = 1, 0o(O) 
finite and = 0 in order to guarantee regularity and 

spatial flatness at the origin, and </>o(°o) = 4 > o , { ( x>) = 0 
in order to ensure asymptotic flatness at infinity as de¬ 
scribed in [lj, (22, |2j, |2J, |23 ; these conditions reduce the 
system (1-3) to an eigenvalue problem for u. The solu¬ 
tion was calculated numerically using finite differencing 
and a shooting routine that searched 1 o. 

The solutions to this set of equations define sequences 
of equilibrium configurations like those shown in Fig. 
1(a). In the curves two important points for each value 
of A are marked: i) the critical point -marked with a 
filled circle- indicating the threshold between the stable 
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FIG. 1: (a) Sequences of equilibrium configurations for dif¬ 

ferent values of A are shown as a function of the central value 
of the scalar field </>o(0). The circles indicate the critical point 
that divides the stable from the unstable solutions. The in¬ 
verted triangles indicate the point at which the binding energy 
is zero. Those configurations between the circles and the tri¬ 
angles (along each sequence) collapse into black holes even 
under infinitesimal perturbations (see (2^] for the full track¬ 
ing of the formation of an event horizon out of an unstable 
boson star). Configurations to the right of the triangles dis¬ 
perse away, (b) The compactness of each solution is shown. 
The critical point is also marked, since the interesting config¬ 
urations would be those boson stars that are stable (to the 
left of the circles). In both plots, the transparent circles indi¬ 
cate the special set of configurations we calculate the emission 
spectrum for in what follows; in the case of A = 0 such con¬ 
figuration is very close to the critical point and it is difficult 
to look at, but it is in there. The configuration marked with a 
square corresponds to a configuration with A = 0, M = 0.499 
and compactness 0.06332; this compactness is the same as 
that of configuration with A = 20 and M = 0.633, see text 
below to take advantage of this remark. 


and unstable branches of each sequence, that is, config¬ 
urations to the left of this point are stable and those to 
the right are unstable as found through the analysis of 
perturbations mm, catastrophe theory 26] and full 
non-linear evolution of the evolution equations [jEl E3 
and ii) the point at which the binding energy Eb = 
M — Nm = 0 marked with an inverted filled triangle, 


where N = f j°d 3 x = J — (j)d v (j)*]d 3 x is 

the number of particles and M = (1 — l/a 2 )r/2 evalu¬ 
ated at the outermost point of the numerical domain is 
the Schwarzschild mass; the configurations between the 
instability threshold and the zero binding energy point 
collapse into black holes whereas those to the right dis¬ 
perse away as shown in |27| . The units for the M and N 
are given in M^Jrn. where M p i is the Planck mass and 
to is the mass of the boson. 

Since the mass of the configurations in Fig. 1(a) 
scales with to, the original use of the self-interaction A 
was to allow bigg er masses even if the mass parameter 
to was fixed [28( and thus BS configurations seemed 
to be -gravitationally- similar to compact objects like 
neutron stars |28| . In the present approach neither m 
nor A will be considered related to any type of possibly 
existing particle, instead they are thought of as free 
parameters that permit faking compact objects. We 
prefer to focus on another property of BSs: the com¬ 
pactness. In Fig. 1(b) the compactness of equilibrium 
configurations is shown. Provided BSs have no defined 
surface we consider that the radius containing 99% 
(i?gg) of the total particle number (following (2§. HU 
where 95% was considered instead) is a reasonable place 
where to measure the gravitational field of the star; 
thus the compactness plotted in Fig. 1(b) is defined 
as 2 A f 9 g/i? 99 , where Ngg is the number of particles 
integrated up to Rgg. One could feel tempted to believe 
that the more compact a BS is more similar to a BH 
seems, however as shown below this happens only in 
terms of how deep the trajectories of test particles can be. 

Time-like geodesics. Given the line element ds 2 = 
— a{r) 2 dt 2 + a{r) 2 dr 2 + r 2 dfi 2 , the equation for time¬ 
like geodesics followed by test particles on the equatorial 
plane reads: 
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( 4 ) 


where L 2 = r 4 tp 2 and E 2 = — a 2 t 2 are the squared angu¬ 
lar momentum and energy at spatial infinity, which are 
the conserved quantities of the test particle related to 
the independence on the azimuthal angle ip and t of the 
space-time respectively; an overdot indicates derivative 
with respect to the proper time of the test particle. The 
geodesics for a Schwarzschild BH are given by 0 with 
the values a 2 = a~ 2 = (1 — Since equation 0 for 

a is linear, we rescale this function so that at infinity 
a(r — > oo ) = 1 /a{r — > oo), thus at infinity a 2 = a~ 2 , 
which implies the coefficient of E 2 in 0 is one for BSs, 
which also happens for the BH. 

At this point it is convenient to draw the effective po¬ 
tential V 2 ff = 4, (l + j for a given angular momen¬ 
tum of the test particle. In Fig. Owe show V 2 ^ for dif¬ 
ferent values of L and compare with the potentials due 
to a BH. It happens that the effective potentials for a BS 
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depend on the compactness of the star only and are in¬ 
dependent of the self-interaction; instead of showing this 
fact we present the potentials for configurations with the 
same mass (M = 0.633 (A/^/to)) and A = 0, 20 with 
compactness 0.1666 and 0.06332 respectively. The loca¬ 
tion of the typical potential barrier for test particles on 
a BS space-time 7| |29j changes from one case to the 
other: the less compact the star (A = 20) the farer the 
location of the barrier from the center. 
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FIG. 2: Here it is shown the effective potential that a particle 
with L 2 = 12 M 2 “feels” for two configurations of the same 
mass M = 0.633 and A = 0, 20, and for a black hole of the 
same mass. At r = 6 M it is possible to see the ISCO for 
a test particle on a BH space-time. The interesting result 
is that for the more compact stars test particles can travel 
deeper inside the star. What is possible to infer is that freely 
falling observers cannot distinguish whether there is a BH or 
a BSs from r = 30M on. 
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FIG. 3: a) The squared denominator of E and L in 0; the 
star chosen is the one with A = 0, central field 0) = 0.271 
and mass M = 0.633(A7/;/m). b) Physical quantities used 
for the calculation of D(r) (see below); the solid line indicates 
these quantities for a black hole of the same mass as the star. 
The results in this plot coincide with those obtained in ||. 


It is known that BSs admit stable circular orbits for 
every radius [§]. Thus, here we do not touch that 
problem and restrict to calculate E and L in terms 
of the metric functions. In order to do so we rewrite 
0) for r = 0 and demand dV/dr = 0 with V = 
[(l + L 2 /r 2 ) — E 2 /a 2 ] /a 2 . This condition implies 


E = 


V a 2 — 


& L = 


(5) 


where a 2 — raa' > 0 all the way as shown in Fig. 0a) for 
a typical BS. Before proceeding to the calculation of the 
emission spectrum we need the angular velocity of a test 

particle given by fl = which is a regular function 

for all values of r. In Fig0b) the values of E, L and the 
useful E — f 1L are shown for the M = 0.633, A = 20 BS 
and a BH of the same mass. 

The accretion disk model is that of a geometrically 
thin, optically thick, steady accretion disk. The power 
per unit area generated by such a disk rotating around a 


central object is 


given by 
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where M is the accretion mass rate and r t is the inner 
edge of the disk, keeping in mind that r is in units of 1 /to. 
For black holes this radius is assumed to be at the ISCO 
(r = 6 M) of the hole. Since BSs allow circular orbits in 
the whole spatial domain we consider that Vi = 0 for BSs. 
Furthermore, assuming it is possible to define a local 
temperature we use the Steffan-Boltzmann law so that 
D(r ) = crT 4 , where a = 5.67 x 10" 5 erg s _1 cm~ 2 I\~ 4 
is the Steffan-Boltzmann constant. Now, considering the 
disk emits as a black body, we use the dependence of T on 
the radial coordinate and therefore the luminosity L(v) 
of the disk and the flux F(y) through the expression for 
the black body spectral distribution: 


L(u) = 4:Trd 2 F(i ') 
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where d is the distance to the source, r % and r/ indicate 
the location of the inner and outer edges of the disk, 
h = 6.6256 x 10 -27 erg s is the Planck constant and 
k = 1.3805 x 10 ~ 16 erg K~ x is the Boltzmann constant 
and i is the disk inclination. The algorithm to construct 
the emission spectrum for such a model of accretion disk 
around a BS and a BH is as follows: 

1) Define the space-time functions a and a by choosing 
one of the equilibrium configurations in Fig. (Tj and 
calculate M. 

2) Define the metric of the equivalent BH through 
otBH = \/l — 2 M/r and cibh = 1/i/l — 2 M/r. 

3) Calculate the angular velocity, angular momentum 
and energy of a test particle for both space-times 
^ BS,BH, LbS,BH > E B S,BH- 

4) Use such quantities to calculate the power emitted in 
both cases Dbs{i") and Dbh(c) defined in (JBJ). 

5) Calculate the temperature of the disk in both cases 

TBs(r) = (Dbs^/ct) 1 / 4 and T B h{t) = (r)/a) 1 / 4 . 

6 ) Use such temperature to integrate the luminosity 0 
Lbs(v) and Lbh(v) for several values of v. 

We choose a configuration with M = 0.633(M 2 ; /m) 
because it has been used in the past (e.g. 0 ), but it could 
be any other; the usual reason to use this configuration 
is that for A = 0 it corresponds to the most compact 
configuration in the stable branch (see Fig.0. We take 
the physical parameters used in [ 8 |, where m = 1.2 x 
10 _ 25 GeV, M = 0.633(M 2 ;/m) = 2.8 x 1O 6 M 0 , i = 60°, 
which we use here too. We used r/ = 30 and 50M and 
found no qualitative difference in the spectrum. As a test 
in the present calculations we show here the results found 
in jl| for such system in Fig. 0 as a particular case, which 
is interpreted as a possible signature of BSs because there 
is a hardening of the spectrum at high frequencies. 

Nevertheless, what we do here is to repeat our algo¬ 
rithm for the same mass of the object but different values 
of A keeping the same value of the mass (the configura¬ 
tions marked with simple circles in Fig. 1). The result 
can be found in Fig. 0] where the spectrum is softened 
when increasing A, in fact we notice that for A = 20 
the spectrum could perfectly be that of the equivalent 
black hole; we dub this particular configuration the Star 
Y and remark its compactness 2 M 99 /.R 99 = 0.06332. 
Under the present model of accretion disk it would be 
rather difficult to distinguish between a black hole and 
a boson star with mass M = 2.8 x 1O 6 M 0 under the 
parameters m = 1.2 x 10 _25 GeV and A = 20 compared 
with the case A = 0 where vL{v) is orders of magnitude 
bigger for the BS than for the BH around the 10 16 iTz 
window. 

Once the algorithm is settled down we try one last 
experiment: we choose a A = 0 star with the same 
compactness as that of star Y. such configuration 
appears in Fig. 1 marked with a square. The mass of 
such star is M = 0.499(Af 2 j/m); in order to achieve a 
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FIG. 4: The emission spectrum from a central star with mass 
M = 0.633 (Mpi/m) and different values of A. The solid line is 
obtained when a black hole is considered as the central object. 
The curve for A = 0 coincides with the one obtained in i8j 
and shows the hardening of the spectrum at high frequencies. 
The case A = 20 is particularly interesting because it seems 
to match the spectrum obtained with the BH. 

mass equal to that of configurations of Fig. 0] we need to 
use m = 0.499/0.633 x 1.2 x 10 _25 GeV. The resulting 
spectrum is shown in Fig. 0 where once again it is 
possible to nearly match the spectrum corresponding to 
that of a BH with the same mass, but this time without 
the need of self-interaction. 
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FIG. 5: The emission spectrum from a disk with a central 
star with mass M = 0.499(M/ ; /m) and A = 0. In order to 
keep M equal to that of configurations in Fig. Hit is necessary 
to set m = 0.499/0.633 x 1.2 x 10 _25 GeV, otherwise we could 
construct a correct spectrum with the incorrect gravitational 
mass. The similarity between the spectrum due to the BS 
and to the BH is not as impressive as that of the Star Y. 
but the hardening of the spectrum is not as strong as that of 
A = 0 in the previous plot. We simply used a star with no 
self-interaction and compactness 2Nqq /Agg = 0.06332 as that 
of Star Y. 

Discussion and conclusions. There are two main 
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results: 1) BSs can play the role of BHCs and 2) the 
emission spectrum is not helpful at determining whether 
or not the BHC is a BS or a BH and in consequence it 
is not enough to say whether the compact object has a 
horizon or not within the accretion model of the present 
approach. There are two ways of faking the spectrum of 
the disk around a BH: i) by using an appropriate value 
of A for a given m and ii) by using an appropriate value 
of to with A = 0. 

We expect the results presented here are disk model 
dependent, however the study of other different accretion 
disk models could show a similar behavior. Nevertheless, 
at this point one can conclude that the spectrum is not 
enough to determine whether a compact object is a BS 
or not, and in any case m and A remain free parameters 
that could help at matching an emission spectrum. 

For instance, in ;9j the accretion of ordinary nucleonic 
gas on black hole candidates is analyzed by assuming 
that the BHC is either a perfect fluid star or a boson 
star instead of a black hole. The result indicates that 
for such cases Type I X-ray bursts should be observed 
because the accreting gas particles pile up on the surface 
of the stars, since in both type of stars the Schwarzschild 
radius lies inside the star. Because such bursts are 
not observed, these two types of stars are ruled out as 


BHCs and it is considered that the BHC is a black hole. 
Nonetheless, the assumption that BSs have surface is not 
correct and a test particle would not be subject to me¬ 
chanical forces due to the presence of a surface; thus the 
analysis in J9j should be redone under such consideration. 

Finally, it could be that under different accretion disk 
models, particular observable BS signatures could appear 
in the electromagnetic spectrum. At the same time, it is 
remarkable the effort invested in the prediction of BS sig¬ 
natures through gravitational wave signals due to parti¬ 
cles surrounding a BS [E[ and due to full 3D perturbations 
@ 
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